Lubricating Bacteria Model for Branching growth of 

Bacterial Colonies 



Yonathan Kozlovsky, Inon Cohen, Ido Golding and Eshel Ben- Jacob 
School of Physics and Astronomy, 
Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69 978, Israel 

Abstract 

Various bacterial strains (e.g. strains belonging to the genera Bacil- 
lus, Paenibacillus, Serratia and Salmonella) exhibit colonial branch- 
ing patterns during growth on poor semi-solid substrates. These 
patterns reflect the bacterial cooperative self-organization. Central 
part of the cooperation is the collective formation of lubricant on 
top of the agar which enables the bacteria to swim. Hence it pro- 
vides the colony means to advance towards the food. One method of 
modeling the colonial development is via coupled reaction-diffusion 
equations which describe the time evolution of the bacterial density 
and the concentrations of the relevant chemical fields. This idea has 
been pursued by a number of groups. Here we present an additional 
model which specifically includes an evolution equation for the lu- 
bricant excreted by the bacteria. We show that when the diffusion 
of the fluid is governed by nonlinear diffusion coefficient branching 
patterns evolves. We study the effect of the rates of emission and 
decomposition of the lubricant fluid on the observed patterns. The 
results are compared with experimental observations. We also in- 
clude fields of chemotactic agents and food chemotaxis and conclude 
that these features are needed in order to explain the observations. 



1 



I. INTRODUCTION 



It is now understood that the study of cooperative self-organization of bac- 
terial colonies is an exciting new multidisciplinary field of research, necessitating 
the merger of biological information with the physics of non-equilibrium processes 
and the mathematics of non-linear dynamics. At this stage, several experimental 
systems have been identified, and preliminary modeling efforts are making sig- 
nificant progress in providing a framework for the understanding of experimental 
observations HTHT^]- 

In nature bacterial colonies must often cope with hostile environmental con- 
ditions. To do so bacteria have developed sophisticated cooperative behavior 
and intricate communication capabilities |14|-- |18|1 . These include: direct cell-cell 
physical interactions via extra-membrane polymers [fL9|,2"0[|, collective production 



of extracellular "wetting" fluid for movement on hard surfaces ]T.21] , long range 



chemical signaling, such as quorum sensing p2|-p4| and chemotactic signalingp] 



25| -E7 , collective activation and deactivation of genes [p8p9| ,3 and even ex- 



change of genetic material p(i|-|32||. Utilizing these capabilities, bacterial colonies 
develop complex spatio-temporal patterns in response to adverse growth condi- 
tions. 

For researchers in the pattern formation field, the above communication mech- 
anisms open a new class of tantalizing complex models exhibiting a much richer 
spectrum of patterns than the models of non-living systems. 

Fujikawa and Matsushita [f|,[3^,[34|] reported for the first time f] that bacterial 



colonies could grow elaborate branching patterns of the type known from the 



1 Chemotaxis is a bias of movement according to the gradient of a chemical agent. 
Chemotactic signaling is a chemotactic response to an agent emitted by the bacteria. 

2 We refer to the first time that branching growth was studied as such. Observations 



of branching colonies occurred long ago [35.36 



study of fractal formation in the process of diffusion-limited-aggregation (DLA) 
f37H39| . This work was done with Bacillus subtilis, but was subsequently ex- 
tended to other bacterial species such as Serratia marcescens and Salmonella 
anatum [ 50 . 



Motivated by these observations, Ben- Jacob et al. conducted new experiments 
with a new species of bacteria that has been isolated from cultures of Bacillus 
subtilis PJ41|J^| . The new species was designated Paenibacillus dendritiformis 



var. dendron [42]. This species is motile on the hard surface and its colonies 
exhibit branching patterns (Fig. 1). The new mode of tip-splitting growth was 
found to be inheritable and transferable by a single cell, hence it is referred to 



as a distinctive morphotype Ji3] , and, to indicate the tip-splitting character of 



the growth, it was denoted T morphotype. In the next section we describe in 
some detail the observations of Ben- Jacob et al. . Additional studies of branching 
colonial growth are reported by Matsuyama et al. and Mendelson and Salhi 



How should one approach the modeling of the complex bacterial patterning? 
With present computational power it is natural to use computer models as a 
main tool in the study of complex systems. However, one must be careful not 



to be trapped in the "reminiscence syndrome", described by J. D. Cowan [45 . 
as the tendency to devise a set of rules which will mimic some aspect of the 
observed phenomena and then, to quote J. D. Cowan "They say: 'Look, isn't 
this reminiscent of a biological or physical phenomenon!' They jump in right 
away as if it's a decent model for the phenomenon, and usually of course it's 
just got some accidental features that make it look like something." Yet the 
reminiscence modeling approach has some indirect value. True, doing so does 
not reveal (directly) the biological functions and behavior. However, it does 
reflect understanding of geometrical and temporal features of the patterns, which 
indirectly might help in revealing the underlying biological principles. Another 
extreme is the "realistic modeling" approach, where one constructs an algorithm 



that includes in details all the known biological facts about the system. Such an 
approach sets a trajectory of ever including more and more details (vs. generalized 
features). The model keeps evolving to include so many details that it loses any 
predictive power. 

Here we try to promote another approach - the "generic modeling" one 
46|, 8 ; 47 , jl| . We seek to elicit, from the experimental observations and the bio- 



logical knowledge, the generic features and basic principles needed to explain the 
biological behavior and to include these features in the model. We will demon- 
strate that such modeling, with close comparison to experimental observations, 
can be used as a research tool to reveal new understanding of the biological 
systems. 

Generic modeling is not about using sophisticated, as it may, mathematical 
description to dress pre-existing understanding of complex biological behavior. 
Rather, it means a cooperative approach, using existing biological knowledge to- 
gether with mathematical tools and synergetic point of view for complex systems 
to reach a new understanding (which is reflected in the constructed model) of the 
observed complex phenomena. 

The generic models can yet be grouped into two main categories: 1. Discrete 
models such as the communicating walkers models of Ben- Jacob et al. |8|,48,|IT 



and the bions model of Kessler and Levine |lg,[i9|]. In this approach, the mi- 
croorganisms (bacteria in the first model and amoebae in second) are represented 
by discrete, moving entities (walkers and bions, respectively) which can consume 
nutrients, reproduce, perform random or biased movement, and produce or re- 
spond to chemicals. The time evolution of the chemicals is described by reaction- 
diffusion equations. 2. Continuous or reaction-diffusion models |50| , |5i| . In these 
models the microorganisms are represented via their 2D density, and a reaction- 
diffusion equation of this density describes their time evolution. This equation is 
coupled to the other reaction-diffusion equations for the chemical fields. In the 
context of branching growth, this idea has been pursued recently by Matsushita 



et al. p2| , Kawasaki et al. [[kJ and Kitsunezaki [S3|. A summary and critique of 
this approach is provided by Rafols |55[| . 

One of the important features of of the bacterial colonies is the lubricant layer 
in which the bacteria swim. A model for the colony should include this feature, 
directly or indirectly |17],5(|. Here we present a model which specifically includes 
the lubricant excreted by the bacteria. The model follows the second approach of 
generic modeling. We represent the various entities: the bacteria, the chemicals 
and the lubricant by continuous fields. 



II. EXPERIMENTAL AND BIOLOGICAL BACKGROUND 

In figure la we show branching patterns of bacterial colonies. Each of these 
colonies is made up of about 10 10 bacteria of the type Paenibacillus dendritiformis 
var. dendron (see for first reference in the literature and |32| for identifica- 

tion). Each colony is grown in a standard petri-dish (8cm in diameter) on a thin 
layer of agar (semi-solid jelly). Figure 2a shows that the branches of the colonies 
have well defined boundary, and the bacteria are confined by this boundary. Fig- 
ures 2b and 2c highlight the constituents of the branches. Figure 2b shows that 
each branch is a layer of fluid on the surface of the agar. Figure 2c shows the 
bacteria, all of which are confined within this fluid. The bacteria cannot move 
on the dry surface and cooperatively they produce a layer of lubrication fluid in 
which they swim. 

Bacterial swimming is a random- walk-like movement, in which the bacteria 
propel themselves in nearly straight runs separated by brief tumbling. Swimming 
can be done only in sufficiently inviscid fluid. To produce such fluid the bac- 
teria secrete lubricant (wetting agents). Other bacterial species produce known 
extracellular lubricants (also known as surfactants, see [0,^,0 and references 
therein). These are various materials (various cyclic lipopeptides were identified) 
which draw water from the agar. The composition and properties of the lubricant 
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of P. dendritiformis var. dendron is not known, but we will assume that higher 
concentration of lubricant is needed to extract water from a dryer agar, and that 
the lubricant is slowly absorbed into the agar (or decomposes). 

In order to move, reproduce and perform other metabolic activities, the bacte- 
ria consume nutrients from the media, nutrients which are given in limited supply. 
The growth of a colony is limited by the diffusion of nutrients towards the colony 
- the bacterial reproduction rate that determines the growth rate of the colony 
is limited by the level of nutrients available for the cells. If nutrient is deficient 
for a long enough period of time, the bacteria may enter a pre-spore state, i.e. 
begin the process of sporulation. They stop normal activity - like movement - 
and use all their internal reserves to metamorphose from an active volatile cell to 
a spore - sedentary durable "seed" . The sporulating bacteria emit a wide range 
of materials, some of which unique to the sporulating bacteria. These emitted 
chemicals might be used by other bacteria as a signal carrying information about 
the conditions at the location of the pre-spores. 

The patterns in figure la are arranged in a diagram according to two control 
parameters: the initial concentration of nutrients (horizontal axis, increasing from 
left to right), and the concentration of the agar, or the dryness of the media 
(vertical axis, increasing from bottom to top). The chirality of the colonies at 
the top row is due to interaction between repulsive chemotaxis (see Sec. |V|) 
and the process of tumbling (see Refs. |48||L6[] for details). It will be ignored in 
this paper. For high level of initial nutrients concentration (right column), the 
patterns are compact, with wide branches. For intermediate levels of nutrients, 
the lower the initial concentration is, the more ramified and less ordered the 
patterns are. The patterns become fractal-like, with fractal dimension decreasing 
for lower levels of nutrients. For the same nutrient level, higher agar concentration 
makes the pattern less dense, with larger gaps between the branches. All of the 
above phenomena could be expected from our knowledge on patterning in non- 
living systems [ 59[- |62 , 16] . Unlike what could have been expected, at the lowest 



concentration of nutrient (leftmost column), the pattern are more ordered with a 
well defined circular envelope. This phenomenon demonstrates the complexity of 
the biological system, and its explanation needs an additional biological feature - 
chemotaxis signaling (see section [Vp. Figure 3 demonstrates that in spite of this 
complexity and the inherit noise in the system, the experiments are controlled 
enough for the patterns to be reproducible. 



III. CONSTRUCTION OF THE MODEL 

The model includes four coupled fields projected on 2D. One field describes the 
density of motile bacteria b(x,t), the second describes the height of the lubricant 
layer in which the bacteria swim l(x,t), third field describes the concentration of 
nutrients n(x,t) and the fourth field is the density of stationary bacteria s(x,t), 
bacteria that enter the pre-spore state. 

We first describe the dynamics of the bacteria and of the nutrient. The two 
react ion- diffusion equations governing those fields have the following form: 

db 

— = movement + reproduction — sporulation (1) 
at 

On 

— = diffusion — consumption 
at 

The sporulation term refers to the transition of motile bacteria into the station- 
ary state, i.e. the pre-spore state. The nutrient diffusion is a simple diffusion 
process with a constant diffusion coefficient. The rate of nutrient consumption is 
proportional to the rate of bacterial reproduction. For the reproduction term in 
([!]) we take: 

reproduction = knb (2) 
where k is a constant rate. This is the usual form included in reaction-diffusion 



models, but see section VI for more detailed discussion. The exact form of the 



sporulation term is not known. For simplicity we take the form: 



sporulation = fib 



(3) 



where the rate \x is constant. 

We now turn to the bacterial movement. In a uniform layer of liquid, bacterial 
swimming is a random walk with variable step length and can be approximated by 
diffusion. The layer of lubricant is not uniform, and its height affects the bacterial 
movement. An increase in the amount of lubricant decreases the friction between 
the bacteria and the agar surface. The term 'friction' is used here in a very loose 
manner to represent the total effect of any force or process that slows down the 
bacteria. It might include, for example, the drag which acts on a body moving 
in shallow layer of viscous fluid. It might include the probability that a flagellum 
will adhere or get tangled with the polymers of the agar. We suggest that the 
bacterial movement depends on the local lubricant height through a power law 
with the exponent 7 > 0: 

movement = V • {D b {l/l M yVb) (4) 

where D b is a constant with dimensions of a diffusion coefficient and Im is the 
height preferred by the bacteria (see below). D b is related to the fluid's viscosity 
and the dryness of the agar might affect this viscosity. But as we argue below, 
this is probably not the main influence of the agar dryness. 
Gathering the various terms gives the partial model: 

^ = V • (D b (l/l M yS7b) + kbn - fib 

—— = D n V 2 n — akbn (5) 
dt w 

ds 

where a is a conversion factor, being the amount of nutrient consumed for repro- 
duction of a new bacterium. The third equation in (|) describes the stationary 
bacteria. Since they are immotile their dynamics include only a source term: 
their conversion from the motile state. 
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The dynamics of the lubricating fluid are also governed by a reaction diffu- 
sion equation. There are two reaction terms: production by the bacteria and 
absorption into the agar. The dynamics of the field are given by: 

?t = -Vj l + f l (b,n,l)-\l (6) 

where Ji is the fluid flux (to be discussed), fi(b, n, I) is the fluid production term 
and A is the absorption rate of the fluid into the agar. 

We assume that the fluid production depends on the bacterial density. As 
the production of lubricant probably demands substantial metabolic efforts, it 
should also depend on the nutrient's level. We take the simplest case where the 
production depends linearly on the concentrations of both the bacteria and the 
nutrients. We presume that the bacteria produce lubricant only when it is needed, 
i.e. when its height is below a threshold height, denoted as l M - We therefore take 
the production term to be: 

fi(b,n,l) =rbn(l M -l) (7) 

where T is the production rate. 

We turn to the flux of the lubricating fluid. The lubricating fluid flows by 
diffusion and by convection caused by bacterial motion. A simple description of 
the convection is that as each bacterium moves, it drags along with it the fluid 
surrounding it. 

J, = -A(ZAm)"VZ + M (8) 

where Di is a constant with dimensions of a diffusion coefficient, Jb is the bacterial 
flux and j is the amount of fluid dragged by each bacterium. The diffusion term of 
the fluid depends on the height of the fluid to the power v > (the nonlinearity in 
the diffusion of the lubricant, a very complex fluid, is motivated by hydrodynamics 
of simple fluids). The nonlinearity causes the fluid to have a sharp boundary at 
the front of the colony, as is observed in the bacterial colonies. The equation for 
the lubricant field is: 
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m 

Of 



V • {D^I/ImYVI + jD b (l/l M yVb) + Tbn(l M - I) - XI 



(9) 



The functional form of the terms that we purposed are simple and plausible, but 
they were not derived from basic physical principles. Therefore we cannot derive 
exact relations between the parameters of those terms and the physical properties 
of the agar substrate. However we can propose some relations between the pa- 
rameters and the agar. In the experiments, the agar concentration is controlled. 
Higher agar concentration gives a drier and more solid substrate. The lubricating 
fluid is composed of lubricant and water. Increasing the agar concentration has 
several effects on the fluid and its dynamics. Since the agar will be drier the ab- 
sorption rate A will increase. It will also make it more difficult for the lubricant 
to extract water from the agar substrate. Therefore the production rate T should 
decrease. More lubricant will be required to extract water so the composition of 
the fluid will change. The lubricant concentration will increase, making the fluid 
more viscous, thus slowing its flow. Therefore the diffusion coefficient Di should 
decrease. A viscous fluid will also slow the bacterial movement, so the bacterial 
diffusion coefficient should also decrease. 

Equation (||) together with equations @ form our model. Before further 
studies of the model we reduce the number of parameters by using dimensionless 
units. We define the new variables: 




(10) 



b' 



bka/ fi, 



n' = nk/n, I' = I/Im 



With the same units we define: 



D' b = D b /D n , D[ = Di/D n , 

V = Tfi/k 2 a, X' = X/fi, f = jn/kalM 



(11) 



Using these variables in (|) and @ and omitting the primes we get: 



8b 
di 



V • (D 6 PV6) +bn-b 
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V 2 n 



bn 



dn 

~dt 
81 

— = V • {D^Vl + jD b FVb) + Tbn{\ -I) -XI 
ds 



(12) 



IV. RESULTS OF NUMERIC SIMULATION 

Figures 4-10 show results of numerical simulations of the model. The figures 
display the sum of the active and stationary bacterial densities, b + s. The 
simulation were done with an explicit method. To reduce the implicit lattice 
anisotropy, a quenched noise was introduced into the diffusion operators. For the 
initial conditions, we set n to have uniform distribution of level no, b to be zero 
everywhere but in the center, and the other fields to be zero everywhere. 

In figure |]we show the fields s, n, b and I. The nutrient field n was consumed 
by bacteria down to a level of n ~ 1 (fi/k in the original units) in the area covered 
by the colony (the nutrient was not completely depleted due to the functional form 
of the bacterial sporulation term, the same as in H| and unlike ]>j|0l)- As the 



nutrient diffuses faster than the other fields it also decreased in the area between 
the colony branches. The field of the motile bacteria b and the lubricant field I are 
overlapping. The motile bacteria are confined to the area covered by lubricant. 
The fronts of both fields have compact support as figure |5] shows. 

In experiments of the bacterial colonies there are two control parameters: the 
agar concentration and the initial nutrient concentration. First we examined the 
effect of changing the latter. As figure |6] shows, we obtained a dense circular 
colony when no was large, a branched pattern when we decreased uq and a DLA- 
like pattern when n was close to 1. Similar effects of decreasing the initial 
nutrient level appear in other reaction-diffusion models |52|-[54],[5(| . 

Changing the agar concentration affects the dynamics of the lubricant fluid. 
We suggested in the previous section that a higher agar concentration relates to a 
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larger absorption rate A and to lower production rate T and diffusion coefficients 
Di and Dj,. It is not a-priori clear what is the exact dependence of each parameter 
on the agar concentration. We use the model to investigate this question. In figure 
|7| we show patterns obtained with different values of the parameters V and A. As 
we expected, increasing A or decreasing T produced a more ramified pattern, 
similar to the effect of higher agar concentrations on the patterns of bacterial 
colonies (figure la). The value of the diffusion coefficient of the lubricating fluid 
Di has less influence on the colony pattern, as can be seen in figure In contrast, 
decreasing the bacterial diffusion coefficient D b produces a ramified pattern, as 
figure | shows. 

In most of the figures we took 7 = 1 and v — 1. There is no a-priori reason 



to take these values, and in figure |10| we show the effect of other values on the 
growth. While the patterns are different, we found that these changes can be 
compensated by other parameters and varying 7 and v has no qualitative effect 
on the conclusions. 

V. CHEMOTAXIS 

The model so far reproduced most of the features of the experimental results 
displayed in figure la, but does not reproduce the transition to ordered patterns 
at the lowest nutrient concentration. We will now extend the lubricant model to 
test for its success in describing this phenomenon. Ben- Jacob et al. suggested 



that this transition is due to chemotaxis and chemotactic signaling |63|- |65| , p!6 
Chemotaxis means changes in the movement of the cell in response to a gradient 
of certain chemical field pB|-|5S|]. The movement is biased along the gradient 



either in the gradient direction or in the opposite direction. Usually chemotactic 
response means a response to an externally produced field, like in the case of 
chemotaxis towards food. However, the chemotactic response can be also to a 
field produced directly or indirectly by the bacterial cells. We will refer to this 
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chemotactic signaling. 
We incorporate the effect of chemotaxis by introducing a bacterial flux J c hem 
due to chemotaxis. The general form is: 

Jchem = C{l)bx{r)Vr (13) 

where r is the chemotactic chemical. x{ r )^ r is the gradient sensed by the bac- 
terium (with x(r) having the units of 1 over chemical's concentration). x(r) is 
taken to be the "receptor law", i.e. x( r ) = K r /(K r + r) 2 |70[| . K r is a constant 
that determines the concentration range of the chemical for which the chemotaxis 
is effective. ((I) (having the same units as a diffusion coefficient) is the bacterial 
response to the sensed gradient (i.e. the effect on the bacterial movement). In 
our model the bacterial diffusion coefficient is D^l", and the bacterial response to 
chemotaxis is ((I) = Co (A>^)- Co is a constant, positive for attractive chemotaxis 
and negative for repulsive chemotaxis. 

Amplification of diffusive Instability Due to Nutrients Chemotaxis: In non- 
living systems, more ramified patterns (lower fractal dimension) are observed for 
lower growth velocity. Based on growth velocity as function of nutrient level 
and based on growth dynamics, Ben- Jacob et al. || concluded that in the case of 
bacterial colonies there is a need for mechanism that can both increase the growth 
velocity and maintain, or even decrease, the fractal dimension. They suggested 
food chemotaxis to be the required mechanism. It provides an outward drift to 
the cellular movements; thus, it should increase the rate of envelope propagation. 
At the same time, being a response to an external field it should also amplify 
the basic diffusion instability of the nutrient field. Hence, it can support faster 
growth velocity together with a ramified pattern of low fractal dimension. The 
bacterial flux due to nutrient chemotaxis is: 

Jnutrient = CnDbV XrJ> , „ " s 9 Vn (14) 

(K n + n) 2 



In figure [TTJ it is shown that as expected, the inclusion of food chemotaxis led 
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to a considerable increase of the growth velocity without significant change in the 
fractal dimension of the pattern. 

Repulsive chemotactic signaling: We focus now on the formation of the fine 
radial branching patterns at low nutrient levels. From the study of non-living 
systems, it is known that in the same manner that an external diffusion field leads 
to the diffusion instability, an internal diffusion field will stabilize the growth. It 
is natural to assume that some sort of chemotactic agent produces such a field. 
To regulate the organization of the branches, it must be a long-range signal. To 
result in radial branches it must be a repulsive chemical produced by bacteria 
at the inner parts of the colony. The most probable candidates are the bacteria 
entering a pre-spore state, which were referred to as the stationary bacteria s. 

As stated above, bacteria may enter a pre-spore state upon starvation. In this 
process they emit a wide range of waste materials, some of which unique to the 
sporulating bacteria. These emitted chemicals might be used by other bacteria 
SIS db SI gnal carrying information about the conditions at the location of the pre- 



spores. Ben- Jacob et al. [p|, [71| , |64| suggested that such materials are repelling the 
bacteria ('repulsive chemotactic signaling') as if they escape a dangerous location. 

The equation describing the dynamics of the chemorepellent contains terms 
for diffusion, production by pre-spores, decomposition by active bacteria and 
spontaneous decomposition: 

dv 

= D r V 2 r + T r s - Q r br - X r r (15) 

dt 

where D r is the diffusion coefficient of the chemorepellent, T r is the emission rate 
of repellent by pre-spores, fl r is the decomposition rate of the repellent by active 
bacteria, and X r is the rate of self decomposition of the repellent. The bacterial 
flux due to repulsive chemotaxis is: 

K r 

Jrepulsive — Cr-^bl ~\_ r h ■ r^Vf (16) 

(K r + ry 

where ( r < (repulsive chemotaxis). 
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In figure the effect of repulsive chemotactic signaling is shown. In the 
presence of repulsive chemotaxis the patterns have a smooth circular envelope, 
while the branches are thinner and radially oriented. 



VI. ANOTHER LOOK ON MODELING BACTERIAL GROWTH 

The model presented in section (|T|) goes along the lines of existing models 



54] , |52|1 in interpreting the bacterial growth terms (reproduction and sporulation). 
There are, however, some biological consideration that may have been overlooked. 
In section |II| we termed knb as 'reproduction' and /i& as 'sporulation'. But looking 
at them term by term unravel discrepancies. The so called "sporulation" term 
in the equations represents, according to the usual interpretation, a constant 
probability per time unit for a bacterium to enter the pre-spore state. This is 
incompatible with our biological knowledge. A nutritional stress is one of the pre- 
requirements for starting the process of sporulation. As long as there is enough 
food for the bacteria to reproduce, there is no nutritional stress. In the model it 
means that as long as the (local) bacterial density increases (i.e. as long as the 
sum of all the reaction terms is positive) there can be no sporulation. 

This leads to a new interpretation of the reaction terms; the resources that the 
bacteria consume from the nutrient (energy and materials) are utilized for two 
main processes: to sustain life and to reproduce. Therefore the reproduction rate 
should be proportional to the nutrient consumed less the amount required for life 
sustaining activities. We assume that the latter is required by each bacterium at 
a constant rate, independent of the nutrient level or the bacterial density. We 
denote that rate by /i, which previously denoted the sporulation rate. Then we 
have instead of (0), with g(b,n) denoting the nutrient consumption rate: 

reproduction = g(n, b) — fib 

if reproduction > 

sporulation = ^ = max (fib — g(n, b), 0) 

-reproduction if reproduction < 
15 



So far we have not changed much the model. All the terms involved in the 
dynamics remain the same, but we give them new interpretations. Only the 
dynamics of the field s changed. Since the dynamics of the other fields are 
decoupled from the field s, the modified model is effectively identical to the 
previous model. 

As we changed the meaning of the bacterial growth terms, we should recon- 
sider their functional form, focusing our attention on the nutrient consumption 
term, g(n, b). The term g(n, b) = knb was taken to be the limit of growth at low 
nutrient concentration, but it is not a-priori evident. We suggest that g should 
have three regimes corresponding to different limits of the nutrient concentration 
and bacterial density. 

1. There is a maximal growth rate of the bacteria, along with a maximal rate of 
nutrient consumption, even for optimal conditions. We denote the maximal 
consumption rate per bacterium by Q n . 

2. When conditions are not optimal and the nutrient concentration is low, it 
might be the limiting factor. Nutrient consumption is then diffusion-limited 
and its rate is proportional to D n n per bacterium, as D n defines the effective 
area from which a bacterium consumes food (assuming no other bacteria 
interfere with this process). 

3. When the nutrient concentration is low and the bacterial density is high, 
there is competition between the bacteria, and the amount of nutrient avail- 
able for each bacterium is proportional to n/b. 

As can be seen in figure 2c, the density of bacteria can be quite high and so 
the third limiting behavior cannot be ignored. The ratio between the nutrient 
diffusion length and the distance between bacteria b~ x l 2 determines which is the 
appropriate regime. We take the nutrient consumption rate to be the minimum 
of the three expressions: 
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g{b,n) = b min(f2 n , kn, en/b) (19) 

where k = D n /a and e is a constant rate (the min function could have been 
replaced by a smooth function with the same limits, but we will put it in this 



way for clarity). We note that the Michaelis-Menten law |7(J of j^-nb as the 
first two expressions as limits, but does not incorporate the third one. The new 
model in dimensionless form is: 

f| = V- {D b PVb)+g(n,b)-b 
at 

dn ^ 2 , 
_ = V n - g{n,b) 

dl 

— = V • {D t FVl + jD b VVb) + Tbn{\ -I) -XI (20) 

C/ V 

ds 

— = m&x(b - g(n,b),0) 

where g(n, b) = b min(— , n, —j) The same transformation of variables (|T0|) and 



lT|) were employed to obtain the dimensionless model. 



Numerical simulations of the above model show (figure [13]) no qualitative 



differences from the model presented in section |V[ It is not the mathematical 
model which is important in this respect, but our understanding of the biological 
system. This emphasizes the differences between the 'generic modeling approach' 
and the approach driving the "reminiscence syndrome" 



VII. CONCLUSION 

We first briefly reviewed experimental observations of branching patterns in 
bacteria of the species Paenibacillus . Both colonial patterns and optical micro- 
scope observations of the bacteria dynamics were presented. 

Our goal in this manuscript was to test a new reaction-diffusion model which 
includes time evolution equation of a lubricant. From a comparison of the model 
simulation and experimental observations we conclude that when a specific bac- 
terial strain is considered, such comparison is not sufficient to tell us if indeed 
the right biological features are included in the model. 
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For more critique test of the models, additional aspects of the growth (such 
as functional dependence of the colonial growth velocity on growth conditions, 
branches size and width distributions etc.) have to be compared with the model 
predictions. One should also compare the theory with more involved experimental 
tests, such as the effect of imposed anisotropy, competition between neighboring 
colonies, and expression of mutants (emergence of sectors) in expanding colonies. 

Our conclusion from the study of bacterial branching growth is that the min- 
imal features of diffusion, food consumption, reproduction and inactivation are 
not sufficient to explain the complete picture of the observed phenomena. We 
believe that additional mechanisms must be introduced, and propose chemotactic 
signaling as plausible one. 

This work has dealt with a continuous model. Such models are not preferable 
to discrete ones. Each has its advantages and disadvantages. The discrete walkers 
model, for example, enables us to include the valuable feature of internal degrees 
of freedom, but is computationally limited in the number of walkers that can be 
simulated, and thus its scaling to the real problem is somewhat problematic. The 
best strategy is to employ in parallel both the reaction-diffusion and the walkers 
approaches. 
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FIGURES 

FIG. 1. a) A diagram of colonial branching patterns of the bacteria Paenibacillus 
dendritiformis var. dendron. The dot at the center of each colony is the initial inoculum. 
The horizontal axis of the diagram is initial nutrient concentration. From left to right 
(in units of gram/liter) it is 0.1, 0.5, 1.0, 2.0 and 3.0 . The vertical axis is agar 
concentration. From bottom to top (in units of Wgram/liter) it is 1.5, 2 and 2.5 
. As the level of initial nutrient concentration is decreased the patterns become less 
organized with fewer branches. However, at the lowest concentration, of 0.1g/l, the 
pattern is ordered with circular envelope, b) An isolated example of branching pattern 
at Q.hg/l pepton and 1.75% agar concentration. 

FIG. 2. Closer look on branches of a colony, a) x20 magnification shows the 
sharp boundaries of the branches. The width of the boundary is in the order of micron, 
b) Numarsky (polarized light) microscopy shows the hight of the branches and their 
envelope. What is actually seen is the layer of lubrication fluid, not the bacteria, c) 
x 50 magnification shows the bacteria inside a branch. Each bar is a single bacterium. 
There are no bacteria outside the branch. 

FIG. 3. Demonstration of the reproducibility of the colonial patterns. Four 
colonies in the same plate (from four inocula). 

FIG. 4. A representative branching pattern produced by the lubricating bacteria 
model. The image shows the four fields of the model. 

FIG. 5. Profile of the fronts of the bacterial field (solid line) and of the lubricating 
field (dashed line) from the 2D model. The fields propagate to the right. Both fronts 
have compact support. 
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FIG. 6. Effect of varying the initial nutrient concentration uq on colony pattern. 
The minimal value of no to support growth is 1. We include the time t it took for the 
colony to grow. 

a. n = 1.1, t = 38089, DLA-like pattern. 

b. n = 1.3, t = 11410. 

c. n = 1.7, t = 3671, branched pattern. 

d. n = 3,t = 721, dense branches. 

e. no = 5,t = 216, a disk. 

The other parameters are: = D\ = 0.5, T = 0.6, A = 0.1, j = 0.01,7 = l,u = 2. 

FIG. 7. Effect of varying A, the fluid absorption rate, on colony pattern. The 
fluid production rate T is 1 in the upper row, and 0.3 in the lower row. In both 
rows A increases from left to right: A = 0.03 (left), A = 0.1 (center), A = 1 
(right) The patterns become more ramified as A increases. Here and in figures [||, 
H decreasing T also produces a more ramified pattern. The other parameters are: 
D b = D l = 1,7 = v = l,n = 1.5, j = 0.01 

FIG. 8. Effect of varying D[, the fluid diffusion coefficient, on colony pattern. Fig. 
(a) and (d) are the same as Fig. (a) and (d). In both rows D\ decreases from left to 
right: D\ = \ (left), D\ = 0.1 (center), Di = 0.01 (right). The patterns do not change 
significantly. 

FIG. 9. Effect of varying Df,, the bacterial diffusion coefficient, on colony pattern. 
Fig. (a) and (d) are the same as Fig. |7] (a) and (d). In both rows Db decreases from 
left to right: Di = 1 (left), Di = 0.1 (center), D\ = 0.05 (right). The patterns become 
more ramified as D b decreases. 
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FIG. 10. Effect of changing the exponents 7 and v. We also print the growth time 
of the colony t. The parameters of the figures are: 

a. j = l )V = l,t = 1249. 

b. 7 = l,v = 2,t = 2069. 

c. 7 = %v = l,t = 4451. 

d. 7 = 2,u = 2,t = 10066. 

The other parameters are: Df, = D[ = 1,Uq = 1.5, T = 1, A = 0.03, j = 0.01. 

FIG. 11. The effect of food chemotaxis on growth. The four patterns differ in the 
values of Q n , the response to the sensed gradient of the nutrient. 
We also print the growth time of the colony t. 

a. C, n = 0,t = 5101 (no chemotaxis). 

b. (n = 10, t = 4957. 

c. Cn = 30, t = 3843. 

d. Cn = 100, t = 2053. 

The other parameters are: = D[ = 1, 7 = v = 1, no = 1.5, T = 0.3, A = 0.1, j = 0. 

FIG. 12. The effect of repulsive chemotaxis on growth. The four patterns differ in 
the values of Cr, the response to the sensed gradient of the chemorepellent. 

a. Cr = (no chemotaxis). 

b. Cr = -20. 

c. £ r = _50. 

d. ^ = -100. 

The other parameters are: Dj, = Di = 1, 7 = v = 1, no = 1.3, T = 0.3, A = 0.1, j = 0. 



FIG. 13. Pattern produced by the modified model of section VI 
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